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Abstract 



We employ a plaquette basis — generated by coupling the four spins in a 2 x 2 
lattice to a well-defined total angular momentum — for the study of Heisenberg 
ladders with antiferromagnetic coupling. Matrix elements of the Hamiltonian 
in this basis are evaluated using standard techniques in angular-momentum 
(Racah) algebra. We show by exact diagonalization of small (2x4 and 2x6) 
systems that in excess of 90% of the ground-state probability is contained in 
a very small number of basis states. These few basis states can be used to 
define a severely truncated basis which we use to approximate low-lying exact 
eigenstates. We show how, in this low-energy basis, the isotropic spin-1/2 
Heisenberg ladder can be mapped onto an anisotropic spin-1 ladder for which 
the coupling along the rungs is much stronger than the coupling between 
the rungs. The mapping thereby generates two distinct energy scales which 
greatly facilitates understanding the dynamics of the original spin-1/2 ladder. 
Moreover, we use these insights to define an effective low-energy Hamiltonian 
in accordance to the newly developed Contractor REnormalization group 
(CORE) method. We show how a simple range-2 CORE approximation to 
the effective Hamiltonian to be used with our truncated basis reproduces the 
low-energy spectrum of the exact 2x6 theory at the ^ 1% level. 



I. INTRODUCTION 



Interest in the study of ladder compounds has been stimulated by suggestions that these 
deceptively simple materials could exhibit some of the critical behavior believed to be re- 
sponsible for high-temperature (high-Tc) superconductivity Moreover, these ladder 
materials — built one chain at a time — could bridge the transition from the one-dimensional 



(ID) chains, where large simulations are carried out routinely, to the two-dimensional (2D) 
structures that are known to be at the heart of the high-Tc compounds 0]. Although this 
might very well be the case, the road from one to two dimensions has been full of (quantum) 
surprises 0. Indeed, it is now known that some of the properties common to the ID- and 
2D-systems — such as the absence of a spin gap — are only shared by the odd-leg ladders. In 
contrast, even- leg ladders have a finite spin gap which should manifest itself in the form 
of an exponential behavior of the magnetic susceptibility at low temperatures; this (expo- 
nential) activation has been confirmed experimentally |Q. Finally, the excitement about 
the ladder materials has recently been fueled even further by some preliminary reports that 
suggest a superconducting transition in some physical realizations of the ladders at a critical 
temperature of about Tc = 12 K p. 

An arsenal of numerical approaches has been employed to elucidate the physics of the 
ladder materials. This includes exact and Lanczos diagonalization techniques 0, quantum 
Monte-Carlo methods |]^, and a density-matrix-renormalization-group approach While 
these techniques have achieved a high degree of sophistication, the selection of a basis has 
not. Indeed, the overwhelming majority — if not all — of the calculations reported in the 
literature have rely on the traditional "Sz" basis. Undoubtedly, the biggest advantage of the 
Sz basis is its simplicity. Indeed, matrix elements of the Hamiltonian are trivially computed 
in this basis. Yet, it is unlikely that the basis represents the "optimal" choice for the 
study of Heisenberg antiferromagnets with isotropic coupling. 

Recent discussions of the possible utility of alternative bases for the study of the ladder 
materials by Martins, Dagotto, and Riera have stimulated us to examine this issue. The 
initial step in our program was to rediscover the "plaquette" basis []10[ and then to examine 



how it might be used in the study of spin ladders. A particularly appealing feature of 
the plaquette basis — generated from the coupling of the four spins in a 2 x 2 lattice to a 
well-defined total angular momentum — is that the states in the basis are the eigenstates of 
the largest two- leg ladder that can be solved in closed form. By this mere fact, it is clear 
that much of the important physics of the problem has been incorporated into the basis. 
Obviously, an important component of our program is the evaluation of matrix elements of 
the Hamiltonian. It is for this that we rely heavily on the sophisticated apparatus of angular- 
momentum algebra [|ll|,|l^ that has been developed in atomic and nuclear physics Hl^JTE 



over many years. 

Our paper has been organized as follows. In Sec. || we describe the plaquette basis and 
compute all the necessary matrix elements of the Hamiltonian. In Sec. |TT| we study small 
2x4 and 2x6 ladders to illustrate the advantage of our approach. In particular, we define 
a severely truncated basis which includes only the 4 lowest, out of 16, 1-plaquette energy 
levels. We use the newly developed Contractor REnormalization (CORE) method to 
construct a low-energy effective Hamiltonian to be used with the truncated basis and find 
that the exact eigenvalues of the 3-plaquette system are very accurately reproduced. Finally, 
we summarize in Sec. 

II. FORMALISM 

The focus of our paper is the plaquette basis. This basis is generated from the coupling of 
the four spins in a 2 x 2 plaquette to a well-defined total angular momentum. To define the 
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basis and to illustrate most of the techniques it is sufficient to concentrate on the 2x4 (or 
two-plaquette) antiferromagnet with isotropic nearest-neighbor coupling J = 1. Moreover, 
we restrict ourselves to the S = 1/2 case, although the formalism can be applied with few 
modifications to arbitrary spin. Indeed, we exploit this flexibility when, in Sec |T|, we discuss 
spin-1 ladders which are in some sense equivalent to the spin-1/2 ladders discussed in this 
section. The model Hamiltonian is given — adopting open boundary conditions — by 

= ^ S, ■ S, = /7o(l) + H,{2) + 2) , (1) 

where denotes nearest-neighbor sites and 

H^{1) = Si • S2 + Si ■ S3 + S2 ■ S4 + S3 ■ S4 , (2a) 
H,{2) = S5 ■ Se + S5 ■ Sr + Se ■ Sg + Sr ■ Ss , (2b) 
V-(1,2) = S3-S5 + S4-S6. (2c) 

Note that this simple two-leg ladder is formed from two interacting n = 4 spin chains; spins 
on the left (right) chain are labeled with odd(even) numbers. Moreover, Hq represents the 
Hamiltonian of an isolated 2x2 system, while V"(l, 2) is the coupling "potential". The 2x2 
Hamiltonian Hq — a simple "antiferromagnet" [|10| — can be re-written in the following form: 



H,{1) = (Si + S4) ■ (S2 + S3) = Li • L2 = i (j? - - L^) , (3) 
where we have defined link and plaquette angular-momentum variables, respectively, as 

Li = Si + S4 ; L2 = S2 + S3 ; Ji = Li + L2 = Si + S2 + S3 + S4 . (4) 

The physically appealing feature of the plaquette basis is that Hq is diagonal in this basis. 
That is, 

^^o(l) 1/1/2, Jimi) = ^[ji(ji + 1) - /i(/i + 1) - hih + 1)] 1/1/2, Ji^i) . (5) 

It is interesting to note that one needs to couple the spins along the diagonal — which are 
the only spins that do not interact — in order to bring the Hamiltonian into a diagonal form. 
The eigenvalues and eigenvectors of the 2x2 Hamiltonian have been listed in Table |. It is 
also instructive to write the ground-state of i^o, namely, the |/i = /2 = l,ji= mi =0) state, 
in terms of the 5"^ basis. That is, 



1 



1, 1 1, 1, 



vj/o) = -= I Tiit) + I itti) - -I TTii) - -I nti) - -I itit) - ol iitt) • (6) 



2' ' 2' ' ' ' 2' 2' 



This simple finding, namely the strong fragmentation of the ground-state probability in the 
Sz basis — but not in the plaquette basis — represents one of our central results. 

The most challenging part of the calculation in the plaquette basis is the computation 
of the matrix elements of the coupling potential 1^(1, 2); recall that in the Sz basis matrix 
elements of the Hamiltonian can be evaluated by inspection. Although by no means trivial, 



the sort of computations which arise in this new basis are done routinely in atomic [13| and 
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nuclear physics |I4| , [T5| . Indeed, over the years a sophisticated formahsm — known generically 
as Racah algebra — has been developed to tackle these computations |Tl],[l2|. It is on these 
techniques that we rely heavily to compute the matrix elements of V{1, 2). We find 



m2)j[,ilsQj2,jrn\V{l,2)\ilJ^^^^ = -( 
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(7) 



Note that we have defined x = ylx+1 and that the +/— sign in the above expression 
should be adopted whenever k + k^k'^k + ji+ ji+ j2 + j'2 = even/odd. Aside from simple 
phases and numerical factors, the matrix elements of V{1, 2) depend on scalar functions (i.e., 
independent of m) known as Racah coefficients; here we cast the matrix elements in terms 
of the more symmetric 6— j symbols |Tl]|T^. Closed-form expressions are readily available 
for the numerical computation of these (and many other) recoupling coefficients [|17|] . 



III. APPLICATIONS 

In the present section we concentrate on small (2x4 and 2x6) ladders to illustrate some 
of the advantages of the new basis. Our aim is to provide convincing evidence, by means 
of a few simple examples, of the utility of the plaquette basis. Hopefully, the evidence will 
be strong enough to persuade some researchers in the field to exploit this basis in future 
calculations. 



A. Distribution of ground-state strength 



In Fig. |1| we display the ground-state probabilities (|(a|\E'o)P with a an element of the 
basis) for the 2x4 Hamiltonian. This picture emphasizes what was already evident from 
Eq. (^, namely, a strong fragmentation of ground-state strength among the many (70) states 
in the Sz basis. Indeed, the only basis states containing a nontrivial amount (~ 15%) of 
strength are the (staggered) Neel states. In contrast, in the plaquette basis a single state 
carries 85% of the ground-state probability and with two states one can practically account 
for the full probability. Moreover, this picture deteriorates little as one increases the size of 
the system. Indeed, for a 2 x 6 system one basis state — out of 132 — carries in excess of 70% 
of the ground-state probability and with only three states one can account for almost 90% 
of it. In Fig. 1^ we show the corresponding distribution of strength for the first excited state 
of the system. Most of the features observed for the ground state remain valid in this case 
as well. Based on this evidence, we believe that the plaquette basis could prove very useful 
in numerical computations of ladder compounds. 
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B. Contractor renormalization group method 



Given that most of the physics of the ground state (and of the low-lying excited states) is 
contained in a very few numbers of plaquette-basis states, it seems natural to attempt some 
form of truncation of the basis so that larger systems may be more readily simulated in the 
computer. The choice becomes obvious upon glancing at Table ||; one should retain the first 
four states, having energies E = —2 and E = —1 (three-fold degenerate), respectively. Note 
that with this choice the link angular momenta — the high-energy degrees of freedom in the 
theory — have been "frozen" at li = I2 = 1, while the plaquette angular momentum, having 
values of ji = 0, 1, becomes the effective low-energy degree of freedom in the new theory. In 
this way the original theory defined on a 16- dimensional Hilbert space (per plaquette) will 
get mapped into a new effective theory having the same low-energy physics as the original 
theory, but defined on a Hilbert space of one fourth the size. The low-energy effective 
theory will now be constructed using the newly developed Contractor REnormalization 



group (CORE) method |T6 



CORE provides a systematic approach at constructing the new low-energy theory using 
contraction and cluster-expansion techniques. The first step into the implementation of 
CORE is the selection of an "elementary" block and a truncation scheme. In our case the 
block is the 2x2 plaquette and the truncation scheme has been described above; the link 
variables become frozen and the plaquette variable is limited to take the values ji = 0, 1. 
Constructing an effective Hamiltonian with the same low-energy properties as the original 
theory — on a system that contains only one single block {Bi) — is straightforward. We obtain. 



{j[m[\Hc{iiBi)\jimi) 



^Ji(ji + 1) 



Shj[Smim[ ■ (8) 



This low-energy Hamiltonian constitutes the range- 1 term in the cluster expansion and is 
denoted by hi{Bi) = Hcs{Bi). One now proceeds to calculate the range-2 contribution to 
the cluster expansion by computing an effective Hamiltonian on a system that contains two 
connected blocks (-Bi and B2). This is the 2x4 Hamiltonian of Eq. (^. In constructing 
the various contributions to the cluster expansion one must pay particular attention to the 
overlaps between the exact eigenstates of the original Hamiltonian and the low-energy basis. 
Since it is simpler to study these overlaps in Hilbert spaces which reflect the symmetries of 
the Hamiltonian, we work with a low-energy basis of definite total spin; recall that [H, J] = 0, 
where J = S = is the total spin of the system. 

We start with the j = sector. In this sector, there are 14 eigenstates of the exact 2x4 
Hamiltonian, while there are only two low-energy states in the truncated basis, namely, 
\(f)i) = \ji=j2 = 0,j = m = 0) and |02) = \ji=j2 = = ttT' = 0). These two states have a 
nonzero overlap with the exact ground state I-Eq^V Indeed, these are the two states that 
dominate the ground-state probability in Fig. |l|. In particular, this implies that both of 
these basis states will "contract" onto the same eigenstate of the Hamiltonian, i.e., 

lim e-*^|0i) oc 1 4°^) and lim e'^^l^a) oc . (9) 

t — >oo t — >oo 

CORE demands that only one low-energy basis state should contract into the ground state. 
CORE also offers a simple solution to this problem: construct a new truncated basis by 
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performing a similarity transformation on the original one so that each state in the new 
basis (^i) contracts onto a unique eigenstate of the exact Hamiltonian, i.e., 

lim e-*^|ei) oc but lim e-*^!^) oc ^) , (10) 

t— >oo t— »oo 

where is the second lowest eigenvalue of the Hamiltonian in the j = sector. This is all 
that is needed. In this way, the relevant matrix elements of the effective Hamiltonian in the 
j = sector become: 



(11) 



(jij2J =0| |i/eff B^) I \J1J2J = 0) = 

' COS 6*0 — sin 6o\ ( Eq^^ ^ \ f '^'^^ sin Oq 
, sin 9q cos 9qJ \ e[^^ J V - sin 6*0 cos Oq 

where E^°^ = -4.293, e[^^ = -2.500, and = -18.482°. 

The construction of the effective Hamiltonian in the j = 1 sector proceeds in a similar 
fashion. In this sector there are 28 eigenstates and three low-energy basis states (each with 
a three-fold degeneracy). These are: |03) = |ji = 1 j2 = 0, j = 1 m), |04) = |ji = 0j2 = l,J = 
Im), and |05) = \ji=j2 = l,j = lm). Since, in this particular case, only the first two states 
need to be transformed, the similarity transformation can again be parameterized in terms 
of a single angle. That is, 

{j[j'j^l\\H,s{BuB2)\\jij2j = l) = 

/cos^i -sin^i 0\ \ / cos^i sin^i 0\ 

sin6li cos6'i e[^^ -sin6li cos6'i 

V ly V V l) 



(12) 



where Ej^'^ = -3.523, e[^^ = -2.915, E^^^ = -2.590, and 9i = 45°. 

Finally, since (up to a five fold degeneracy) there is a unique j — 2 state in the low energy 
basis, namely, l^g) = \ji—j2 — i,j— 2m), the effective Hamiltonian in this sector is simply 
given by 

{j[j',j^2\\H,s{Bi,B2)\\hj2j^2) = = -2.207 . (13) 

Collecting all the above results we can now write the effective 2x4 low-energy Hamiltonian 
in the direct product basis: 

Yl(jij2j\\Hes{Bi, B2)\\jij2j){j[m[, j'^m'^lj m[ + m'^) {jirrii, j2m2\3 rrii + ma) , (14) 
j 

where (j'lmi, jamaljm) are Clebsch-Gordan coefficients. The range-2 contribution to the 
cluster expansion is obtained by simply removing from ifeff(-Bi, -B2) those range-1 terms 
that have already been included in the single- block calculation, i.e., 

/i2(5i, B2) = H,s{Bi, B2) - hi{B^) - h^{B2) . (15) 

To construct the renormalized Hamiltonian one must continue this procedure, indefi- 
nitely, on larger and larger connected blocks. Here we will stop at the range-2 contribution. 
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Note that for our choice of basis, this range-2 approximation already takes into account 
correlations among next-to-next nearest neighbors. The approximate (up to range-2) renor- 
malized Hamiltonian becomes 



H,en = E + h2iBj, . (16) 



It is instructive to use this approximation to compute the low-energy spectrum of the 2x6 
(three-plaquette) Heisenberg antiferromagnet. For this case, the range-2 approximation 
yields 

Hren = hi{Bi) + hi{B2) + /ll(53) + h^iB^, B2) + h2{B2, B3) . (17) 

This expression is useful as it suggests when the cluster expansion might become rapidly 
convergent. If an optimal basis has been chosen — provided that one exists — one might hope 
that most of the low-energy spectrum could be generated by the range- 1 terms, leaving the 
range-2 terms in charge of the fine tuning. Alternatively, an optimal basis could generate — 
dynamically — two energy scales in the problem; a large one associated with physics within 
the blocks and a small one associated with the "residual interaction" between the blocks. 
Our investigations show that — for our truncated basis — such is indeed the case. First, one 
can simply compare the individual matrix elements of hi and /i2. Those of the former are 
typically 4 to 8 times larger than those of the latter. Second, one can arrive at the same 
conclusion by mapping the original isotropic spin- 1/2 ladder onto an equivalent anisotropic 
spin-1 ladder as follows. Note that for the 2-plaquette system [see Eq. (^] 

Hoil) + Ho{2) = (Si + S4) • (S2 + S3) + (S5 + Ss) ■ (Se + St) = Li ■ L2 + L3 • L4 . (18) 

In the truncated basis, all link angular momenta are 1; hence the Li are spin-1 operators. 
Now consider ^(1,2) = (S3 ■ S5 -|- S4 ■ Sg). In the truncated basis, symmetries permit the 
interaction to be written as 

V(1,2) = S3-S5 + S4-S6 

^ ^[(82 + S3) ■ (Ss + Ss) + (Si + S4) ■ (Se + St)] 

= i[Li-L4 + L2-L3] . (19) 

Hence, in the truncated basis — which carries much of the important physics of the problem — 
the isotropic 2-plaquette spin-1/2 system is equivalent to a 1-plaquette spin-1 system where 
the coupling along the rungs (Li ■ L2 + L3 ■ L4) is four times as strong as the coupling 
between the rungs (|[Li ■ L4 -|- L2 • L3]). In this context our basis is optimal in the sense 
that the much stronger "rung" couplings are diagonal in it. We also note that this mapping, 
in conjunction with CORE techniques similar to those discussed above, is likely to provide 
a useful starting point for formulating a renormalization group transformation which can 
permit us to estimate properties of infinite ladders. This will be the topic of a future 
publication. 

In Table |I| we display the low-energy spectrum of the 2x6 Hamiltonian using a variety 
of approximations; the ratio to the exact value appears in parenthesis. Recall that the 
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spectrum has been computed with open boundary conditions. The states have been classified 
according to their total spin, which is listed in the first column. In the second column we 
report the results from an extreme "weak-coupling" calculation. In this approximation the 
2x2 plaquettes are treated exactly but the residual interaction between the plaquettes is 
neglected [i.e., V{j,j-\-l) = 0; see Eq. (^]. We observe that at the 80 to 90 percent level, the 
spectrum is, indeed, accounted for by the mere selection of the basis. In the third column we 
report a calculation which uses the exact 2x6 Hamiltonian but with the truncated low-energy 
basis. Departures from the exact results are only 5 to 10%. In principle, this truncation 
could enable the simulations of larger systems, as the number of states increases with the 
number of plaquettes [Np) only as 4^^, rather than as However, in many applications, 
this level of accuracy may still be insufficient. Amazingly, a dramatic improvement on these 
calculations results from expending the very little additional effort required to construct 
the CORE effective Hamiltonian for the truncated basis. As shown in Table |I|, CORE- 
improved calculations yield results that range from a fraction of 1% to a few percent of the 
exact answer. Note that these results were obtained by the diagonalization of — at most — 
20 X 20 matrices; instead, the exact calculation in the 5*2 basis requires a diagonalization 
of a 924 X 924 matrix. In the near future, we plan to use this renormalized Hamiltonian, 
perhaps including range-3 contributions, to simulate larger systems. 

We conclude this section with a brief comment about the doping of the ladders. In 
order to gain some qualitative insight into the nature of hole correlations, Dagotto and 
collaborators introduced two energy scales in the problem: a large exchange coupling J' 
along the rungs relative to a small exchange coupling J along the chains [0]. They observed 
that if a pair of holes is added to the system, the energy will get minimized whenever 
the two holes go into the same rung in the ladder. Further, they concluded by means of 
numerical evidence, that most of the arguments developed for the anisotropic case remain 
valid even when J' ^ J; note that to date, the physical realization of the ladders seem to 
obey the isotropic relation In this paper we have only considered the isotropic case. A 
particularly gratifying aspect of the mapping described above is the natural appearance in 
the isotropic system of distinct energy scales in the effective spin-1 ladder. This appears to 
permit carrying over the arguments of Dagotto and collaborators with little modification. 
Now it is the rungs on the spin-1 ladder which are strongly bound, consisting of a pair of 
"frozen" spin triplets coupled to an overall angular momentum of zero; interactions between 
rungs are relatively weak. In terms of the original spin-1/2 ladder, we may conclude that 
individual plaquettes interact only weakly. Hence the weak-coupling limit defined above 
should be a reasonable approximation when applied to very large systems just as we have 
found it to be by comparing with exact results for the 3-plaquette case. In this limit it 
is simple to see that the ground-state of the system consists of all plaquettes being in the 
lowest j = state, with energy per plaquette of —2 (see Table |). When a pair of holes is 
introduced into the system the holes can go into two different plaquettes at a cost in energy 
of +2 (the lowest energy of three spins in a plaquette is equal to —1). Alternatively, the holes 
can go into a rung — or along a chain — in the same plaquette at a cost in energy of +5/4 (the 
lowest energy of two spins in a plaquette is equal to —3/4). Hence, it becomes energetically 
favorable for the two holes to bind and break as few j = plaquettes as possible. Moreover, 
as another pair of holes is added into the system, it becomes energetically favorable — at 
least for holes with no mobility — for the four holes to go into the same plaquette, rather 
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than for the new pair to break another plaquette. 

IV. CONCLUSIONS 

Wc have employed a plaquette basis for the study of antiferromagnetic Heisenberg lad- 
ders. The states in the basis represent the eigenstates of the 2x2 plaquette and are con- 
structed from the angular-momentum coupling of the four spins in the plaquette. Matrix 
elements of the Hamiltonian were computed in this basis and were expressed in terms of 
a product of five Racah coefficients. These expressions are considerably more complicated 
than the corresponding ones obtained using the conventional Sz basis. Yet, they can be 
efficiently computed by employing angular-momentum techniques that have been developed 
over the years in atomic and nuclear physics. Moreover, this basis seems to capture some of 
the important physics of these complicated systems. Indeed, we have shown that the distri- 
bution of ground-state and first-excited-state strength is concentrated in a very few number 
of states. This is in contrast to the Sz basis where the strength is strongly fragmented. This 
concentration of strength among a few states provides a very natural truncation scheme 
for the basis. We selected a low-energy basis which reduces the size of the Hilbert space 
by a factor of four per plaquette. By using the original (2 x 6) Hamiltonian in this trun- 
cated space we obtained a low-energy spectrum that was within 10% of the exact answer. 
However, by improving the Hamiltonian — via CORE — we were able to get within 1% of the 
exact answer. Moreover, the mere selection of the basis dynamically generates two-energy 
scales in the problem — even in the case of an isotropic coupling. This phenomenon is most 
clearly understood by mapping the original isotropic spin-1/2 ladder onto an effective spin-1 
ladder in which the coupling along rungs is four times stronger than the coupling between 
rungs. This separation of scales is important for the development of qualitative insights into 
the nature of hole doping. In particular, it supports the notion that holes will go into the 
ladders in such a way as to disturb the minimum number of j = plaquettes. 

In summary, we have provided solid evidence in support of a plaquette basis for the 
study of Heisenberg ladders. Much work remains to be done, such as the implementation of 
a Lanczos diagonalization procedure and the study of larger systems using the renormalized 
range-2 (or even range-3) CORE Hamiltonian. Yet, we believe that the mere selection of 
a basis could play a prominent role in the elucidation of the important physics behind the 
ladder materials. 
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FIGURES 



FIG. 1. Distribution of strength for the ground state of the 2x4 Hamiltonian in the Sz basis 
and in the plaquette basis. 

FIG. 2. Distribution of strength for the first-excited state of the 2x4 Hamiltonian in the Sz 
basis and in the plaquette basis. 
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TABLES 

TABLE L Eigenvalues and eigenvectors of the 2x2 system. 



^1 I2 ji E degeneracy 

1 i ^2 i 

111-1 3 

0000 1 

0110 3 

1010 3 

112+1 5 



TABLE IL Low-lying spectrum for the 2x6 Heisenberg ladder. The different approximations 
are explained in the text. Quantities in parenthesis represent the ratio to the exact value. 



s 


Weak Coupling 


Truncated 


CORE (Range-2) 


Exact 





-6.000 (0.909) 


-6.335 (0.959) 


-6.588 (0.998) 


-6.603 




-4.000 (0.818) 


-4.396 (0.899) 


-4.899 (1.002) 


-4.888 




-4.000 (0.836) 


-4.250 (0.889) 


-4.882 (1.021) 


-4.783 


1 


-5.000 (0.844) 


-5.580 (0.942) 


-5.928 (1.001) 


-5.924 




-5.000 (0.909) 


-5.112 (0.930) 


-5.458 (0.993) 


-5.498 




-5.000 (0.983) 


-4.666 (0.917) 


-5.132 (1.009) 


-5.087 




-4.000 (0.797) 


-4.538 (0.905) 


-5.092 (1.015) 


-5.017 




-4.000 (0.850) 


-4.125 (0.877) 


-4.809 (1.022) 


-4.705 


2 


-4.000 (0.823) 


-4.413 (0.908) 


-4.862 (1.000) 


-4.862 




-4.000 (0.925) 


-3.875 (0.890) 


-4.426 (1.017) 


-4.352 


3 


-3.000 (0.895) 


-2.750 (0.821) 


-3.414 (1.019) 


-3.351 
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